A Composite Fermion Description of Rotating Bose-Einstein Condensates 
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We study the properties of rotating Bose-Einstein condensates in parabolic traps, with coherence 
length large compared to the system size. In this limit, it has been shown that unusual groundstates 
form which cannot be understood within a conventional many-vortex picture. Using comparisons 
with exact numerical results, we show that these groundstates can be well-described by a model 
of non- interacting "composite fermions". Our work emphasises the similarities between the novel 
states that appear in rotating Bose-Einstein condensates and incompressible fractional quantum 
Hall states. 
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It has proved fruitful in fractional quantum Hall sys- 
tems to account for the many-body correlations in- 
duced by electron-electron interactions by introducing 
non-interacting "composite fermions" Recently a 

similar approach has been employed to show that the cor- 
related states arising from interparticle interactions in di- 
lute rotating confined bose atomic gases can be described 
in terms of the condensation of a type of composite boson 
1^. Here, we demonstrate that a transformation of the 
system of rotating bosons to that of non-interacting com- 
posite fermions is also successful in accounting for these 
correlated states. Our results establish a close connec- 
tion between the groundstates of rotating confined Bose- 
systems and the correlated states of fractional quantum 
Hall systems Q. 

While the trapped atom gases have been shown to 
Bose-condense the response of these condensates 

to rotations has not, as yet, been measured experimen- 
tally. Theoretically, it is clear that there exist various 
different regimes. Within the Gross-Pitaevskii frame- 
work, which requires macroscopic occupation of the sin- 
gle particle states, the system forms vortex arrays at both 
long 1^ and short coherence lengths (compared to 
the size of the trap), which are reminiscent of Helium- 
4. Here, following Ref. we choose to study the sys- 
tem in the limit of large coherence length without de- 
manding macroscopic occupation numbers. This allows 
us to study both the regime considered in Ref. ||^ , as well 
as regimes of higher vortex density where the quantum 
mechanical nature of the vortices will be most preva- 
lent. Indeed, in Ref. |3| it was shown that, in general, 
the groundstates of the rotating boson system cannot 
be described within a conventional many-vortex picture. 
Rather, the system was found to be better described in 
terms of the condensation of "composite bosons" - bound 
states of vortices and atoms - across the whole range of 
vortex density. In the present paper, we show that a de- 
scription in terms of non-interacting composite particles 
with fermionic statistics also provides a highly accurate 
description of the rotating bose system: specifically, it 
enables us to predict many of the features in the en- 
ergy spectrum and to form good overlaps with the exact 



groundstate wavefunctions. In addition, this description 
indicates a close relationship between the properties of 
rotating Bose systems and those of fractional quantum 
Hall systems. 

In a rotating reference frame, the standard Hamilto- 
nian for N weakly interacting atoms in a trap is 
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where we have used the trap energy, hy^ K/m = Tiluq as 
the unit of energy and the extent, {fi^ /MKY/'^, of the 
harmonic oscillator ground state as the unit of length. 
{M is the mass of an atom and K the spring constant 
of the harmonic trap.) The coupling constant is defined 
as ?7 = ATTna{h? /MK)^^/'^ where n is the average atomic 
density and a the scattering length. The angular velocity 
of the trap, oj, is measured in units of the trap frequency. 

Throughout this work, we make use of the limit of weak 
interactions (77 ^ 1). It was shown in Ref. [|j that in this 
limit the system may be described by a two-dimensional 
model with a Hilbert space spanned by the states of 
the lowest Landau level: ipm{r) oc z™ exp(— zz*/2), 
where m is the angular momentum quantum number 
(m ~ 0,1,2...) and z = x + iy. The kinetic energy 
is quenched and the groundstate is determined by a bal- 
ance between the interaction and potential energies. Not- 
ing that the ^-component of the angular momentum, L, 
commutes with the Hamiltonian, the total energy, scaled 
by rj, may be written 



E/ij^VN{L) + il-u;)/f]L, 



(2) 



where Vn{L) is the interaction energy at angular momen- 
tum L. While this separation holds for all energy eigen- 
states, we choose Vn{L) to denote the smallest eigenvalue 
of the interactions at angular momentum L. Since the 
interactions are repulsive, Vn{L) decreases as L increases 
and the particles spread out in space; a tendency that is 
opposed by the term (1 — uj)/riL describing the parabolic 
confinement. Thus, as the rotation frequency lo is varied, 
the groundstate angular momentum will increase, from 
L = at Lj = 0, to diverge as ^ 1 (when the trap 
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confinement is lost); our goal is to describe the sequence 
of states (of different L) through which it passes. 

We have obtained the groundstate interaction energies, 
Vn{L), for = 3 to 10 particles, from exact numerical 
diagonalisations within the space of bosonic wavefunc- 
tions in the lowest Landau level While the interac- 
tion energy Vn{L) does decrease with increasing angular 
momentum, it is not a smooth function of L. Thus, the 
groundstate angular momentum, obtained by minimis- 
ing (^), is not a smoothly increasing function of ui. As 
shown in Fig. certain values of angular momentum, 
corresponding to downward cusps in Vn{L), are partic- 
ularly stable, and are selected as the groundstate over a 
range of u. 
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FIG. 1. Groundstate angular momentum as a function of 
rotation frequency, ui, for A'^ = 3 ^ 10 particles 

The existence of certain angular momentum states of 
enhanced stability is reminiscent of the "magic-values" 
of angular momentum for electrons in quantum dots 
|pO|-p^ ; by analogy we will also refer to the stable angular 
momenta of the bosons as the magic values. Indeed, the 
system of bosons we study is precisely the bosonic vari- 
ant of the (fermionic) problem of a parabolic quantum 
dot in strong magnetic field, with ujq playing the role of 
the magnetic field and {l — uj)/r] the role of the parabolic 
confinement, and with (5-function interactions replacing 
the more usual Coulomb repulsion. As we shall explain, 
the magic values of angular momenta for the bosonic and 
fermionic systems are, in fact, closely related. This is a 
corollary of our principal result, to which we now turn, 
that much of the structure appearing in Fig. |l| can be 
interpreted simply in terms of the formation of bound 
states of bosons and vortices behaving as non-interacting 
composite particles with fermionic statistics - "compos- 
ite fermions"(CF). 

It is known that, for homogeneous systems, interact- 
ing bosons and interacting fermions within the lowest 
Landau level have many features in common. For exam- 



ple, there exist certain filling fractions |Q of both the 
boson and fermion systems at which interactions lead 
to incompressible groundstates, with wavefunctions that 
may be related by a simple statistical transformation if 
l/vp = + 1 1 15 1 (i^B and vp are the filling frac- 

tions of the bosons and fermions). These similarities 
arise from the remarkable effectiveness of mean-field ap- 
proximations to Chern-Simons theories of such systems 
Q . Here, we are interested in an inhomogeneous system, 
in which the bosons are subject to a parabolic confine- 
ment. Jain and co-workers [|l^,^^ have shown that the 
fermionic equivalent of this problem - interacting elec- 
trons in a quantum dot - can be well-described in terms 
of properties of non-interacting composite-fermions. Mo- 
tivated by the successes of their theory, we apply a similar 
transformation to describe the present bosonic problem. 

Specifically, we make the following ansatz for the 
many-boson wavefunction 



^ansatz 



(3) 
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where ^'£^^({zi}) is a wavefunction for some fermionic 
particles - the composite fermions. Multiplication of 
the antisymmetric CF wavefunction by the Jastrow pref- 
actor generates a completely symmetric bosonic wave- 
function. V projects the wavefunction onto the low- 
est Landau level, which amounts to the replacement 

Zi™ {^^r"™ (n > m) ; (n < m)} for all terms in 
the polynomial part of the wavefunction. For a full dis- 
cussion see Ref. |Q . (For ease of presentation, we omit 
exponential factors and normalisation constants from all 
wavefunctions.) 

The transformation (|^) causes the boson wavefunction 
to describe a half vortex around the position of each other 
particle in addition to the motions described by ^f*-^^. 
One can therefore interpret a composite fermion as a 
bound state of a boson with a half vortex (cf . Ref. |^ ) . 
As a result, the angular momentum of the bosons, L, is 
increased with respect to that of the composite-fermions, 
LcF, according to 



L = LcF + N{N -l)/2. 



(4) 



Note that the transformation (||) relates l/vcF = ^/vb — 
1 and is not the same as that used in Ref. |l5). There 
are an unlimited number of fermion <-> boson mappings 
that one can effect through transformations of the form 
(^). In Fig. H we present a schematic of how, by subse- 
quent attachments of half- vortices - each causing an ad- 
dition of N{N —l)/2 to the angular momentum - one can 
transform from the composite bosons (CB) introduced in 
Ref. 1^ to the composite fermions used here (CF), to the 
bare boson system in which we are interested (B), and 
finally to a fermion system (F). 
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FIG. 2. A schematic relating the transformation we are 
performing to those that have been used previously. 

We introduce the fermion system (F) to point out that 
the composite fermions (CF) which we use to describe 
the boson sy stern (B) are the same as those used by Jain 
and Kawamura |13| to describe interacting electrons in 
quantum dots (F). The predictions of the energy spec- 
trum flowing from a model of non-interacting compos- 
ite fermions will therefore be identical in the boson and 
fermion systems up to the shift Lp = L + N{N — l)/2. 

In the spirit of Ref. we shall consider the CFs, de- 
scribed by '^1^^{{zi}), to be non-interacting, and look 
at the variation of the minimum kinetic energy of the 
CFs as a function of the total angular momentum. We 
further assume that a composite fermion in the Landau 
level state {n, m) (with Landau level index n = 0, 1, 2 . . ., 
and angular momentum m = —n, —n + 1 . . .) has an en- 
ergy En = {n + l/2)EcF, where Ecf is some effective 
cyclotron energy. These assumptions may be viewed as 
a mean-field treatment of the appropriate Chern-Simons 
theory for this system; ultimately, they are justified by 
the predictive successes of the resulting theory. 

Figure || shows the resulting groundstate energy of 
non-interacting CFs as a function of L = Lcf + N{N — 
l)/2 for N = 7,8,9, together with the exact interaction 
energies Vn{L). 
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FIG. 3. The circles show the groundstate energies of 
non-interacting composite fermions as a function of angular 
momentum, for A'' = 7, 8, 9 (filled circles identify the sets, L*pf, 
of angular momenta at which the CF energy has a downward 
cusp). The exact interaction energies, Vn{L) are shown as 
solid lines for comparison. (All curves are offset for clarity.) 



It is apparent that the composite fermion energies fail 
to capture the rapid rise in the exact energies at small an- 
gular momenta; this can be interpreted as a failing of the 
assumption of a constant effective cyclotron energy Ecf- 
The principal success of this approach is the identifica- 
tion of the cusps in the exact energy Vn{L): at almost all 
of the angular momenta for which the composite fermion 
kinetic energy shows a downward cusp (we label these 
sets of angular momenta by i^), there is a correspond- 
ing cusp in the exact energy Vn{L). Since a "magic" 
angular momentum of the boson system must coincide 
with a downward cusp in V]\[{L), the set represents a 
set of candidate values for the magic angular momenta. 
For example: for N — 7, the CF model predicts all nine 
actual magic numbers and in addition identifies cusps 
which do not become groundstates for a further three 
L e L^. These missing values are not necessarily a fail- 
ing of the composite fermion model. The main failing, to 
which we will return later, is that there is a small number 
of magic angular momenta that are not identified. 

Not only does the composite fermion model success- 
fully identify the majority of the magic angular momenta, 
as we now show it also provides a very accurate de- 
scription of the associated wavefunctions. The composite 
fermion wavefunctions corresponding to the angular mo- 
menta are the "compact states" discussed in Ref. . 
For these states, the composite fermions occupy the low- 
est available angular momentum states within each Lan- 
dau level. As an illustration, for A^ = 4, there is a cusp 
in the composite fermion energy at L = 8 {Lcf = 2), 
at which the composite fermions occupy the single par- 
ticle states (n,m) = {(0, 0), (0, 1), (0, 2), (1, -1)}. The 
wavefunction VP?^ is formed as a Slater determinant 
of these states, and the bosonic wavefunction \]/|^"''"*^ is 
constructed via Eq.(^. 

In the cases L — and L — N{N — 1), this procedure 
yields the exact groundstate wavefunction for all A^. At 
L = 0, there is only one many-body state within the low- 
est Landau level (all bosons occupy the to = state) ; the 
ansatz (^) has non-zero overlap with this state, so must 
(trivially) be the groundstate. For L — N{N — 1), the 
lowest energy composite fermion state is formed from the 
states {(0,0), (0, 1), ....(0, A^)}. The Slater determinant 
of these states may be written ^f*^^ = ni<j(-^i ~ ^j)^ 
which, inserted in (^), generates the bosonic Laughlin 
state ^r'^"'"^*^ = ni<j(^i - Zj)"^- (This state is in the low- 
est Landau level, and projection is unnecessary.) Since 
this wavefunction vanishes for Zi — Zj (i ^ j), it is the 
exact zero energy eigenstate of the (5-function two-body 
interaction potential. 

At intermediate values of the angular momentum, our 
ansatz (^) is not, in general, exact. We have performed 
numerical calculations to determine the overlaps of the 
ansatz wavefunctions with the exact groundstate wave- 
functions, |^\[/2"''''*^|vl'|^'''=*)|. We list these overlaps in 
Table || at each of the angular momenta, L G L^, se- 
lected by the non-interacting composite fermion model. 
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In general, the ansatz (||) has an overlap of close to 
unity with the exact groundstate: the composite fermion 
model provides an excellent description of these states. 
Small overlaps can occur when the composite fermion 
model does not produce a unique ansatz ~ i.e. when two, 
or more, sets of single particle states for the composite 
fermions have the same kinetic energy at a given L (e.g. 
N = 6, L = 12). In these cases, the overlaps could be 
improved by diagonalising the Hamiltonian within the 
space of states spanned by the two ansatz states. 

Owing to the impressive agreement between the ansatz 
wavefunctions and the exact groundstates at L^, an 
accurate description of the groundstate angular momen- 
tum as a function of rotation frequency can be obtained 
using only this set of ansatz wavefunctions. Minimizing 
the expectation value of the energy (||) within this set of 
ansatz wavefunctions, one obtains a groundstate angular 
momentum as a function of {1 — uj)/ri that is in excellent 
agreement with the exact results shown in Fig. |l|. This 
approach does, however, omit a small number of magic 
angular momenta. In some cases, these are magic values 
identified by the composite fermion model, but for which 
the expectation value of the energy happens not to be 
sufficiently low to become stable (TV = 6, L = 12; TV = 9, 
L = 33,37; N — 10, L — 38). The most important omis- 
sions are the magic angular momenta at7V = 8,L=12, 
N — 10, L = 16, 21 for which there are no features in 
the composite fermion kinetic energy that would suggest 
a stable angular momentum state. We believe that this 
emergent structure at larger numbers of particles repre- 
sents many-body correlations that are not captured by 
the non-interacting composite fermion model used here. 
(Some of these states are correctly identified by the com- 
posite boson approach Q.) They could be related to 
the incompressible states, such as v = 4/5, of quan- 
tum Hall systems which cannot be explained in terms 
of non-interacting composite fermions alone, but require 
an additional 'particle-hole' transformation. This view is 
strengthened by the observation that related magic an- 
gular momenta also appear in the exact groundstate en- 
ergy of electrons in quantum dots interacting by Coulomb 
forces^p to the the shift Lp = L + N{N - l)/2 {e.g. 
Ref. ||l^ identifies a stable state of = 10 electrons at 
Lp = 61 - equivalent N — 10, L = 16 of the present 
bosonic model). The study of this additional structure is 
beyond the scope of the present work. 

In summary, we have studied the properties of rotat- 
ing Bose systems in parabolic traps in the limit of large 
coherence length. Through comparisons with exact re- 
sults for small systems, we showed that many of the fea- 
tures of the exact spectrum of the bosons can be un- 
derstood in terms of non-interacting composite fermions. 
The non-interacting composite fermion model leads to 
(1) the identification of a set of candidate values for the 
stable angular momenta of the bosons, and (2) associ- 
ated many-body wavefunctions that have large overlap 
with the exact groundstate wavefunctions. The successes 
of the mapping to composite fermions indicate that the 



groundstates of rotating Bose-Einstein condensates, in 
the limit of large coherence length, are closely related 
to the correlated states appearing in fractional quantum 
Hall systems. 
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TABLE I. For each number of bosons. A'', the angular mo- 
menta, L%, at which the non- interacting composite fermion 
description predicts a downward cusp are given, together with 
[in brackets] the overlaps of the ansatz wavefunction with 
the exact groundstate wavefunction. Where a given angular 
momentum appears more than once for fixed N, the compos- 
ite fermion model provides more than one candidate ground- 
state. 
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